Wide-angle perfect absorber/thermal emitter in the THz regime 
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We show that a perfect absorber/thermal emitter exhibiting an absorption peak of 99.9% can 
be achieved in metallic nanostructures that can be easily fabricated. The very high absorption is 
maintained for large angles with a minimal shift in the center frequency and can be tuned throughout 
the visible and near-infrared regime by scaling the nanostructure dimensions. The stability of the 
spectral features at high temperatures is tested by simulations using a range of material parameters. 
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Since the beginning of the last century it is known that 
a perfect thermal emitter follows Planck's law of black- 
body radiation[l]. Realistic structures, however, gener- 
ally do not follow Planck's law, but exhibit a smaller 
emission. The properties of these emitters strongly de- 
pend on the materials and their shapes. From the absorp- 
tion spectra of a structure the emission properties can be 
deduced since Kirchhoff 's law directly relates the absorp- 
tion with the emissivity. The emission is then determined 
by multiplying the emissivity with the black-body radi- 
ation spectrum. Using Photonic Crystals [3, y, IJjit has 
been shown that this approach is also valid for periodi- 
cally structured materials. For a number of applications 
such as thermo-photovoltaic converters, it is necessary to 
control the spectral properties to achieve, e.g., selective 
emitters in a narrow frequency band corresponding to 
the band gap of solar cells [5|,|6|]. In the case of structured 
metallic surfaces, the changes in the emission spectra are 
based on surface waves coupled to the external radia- 
tion through the periodic surface. 0, H, Q Alternatively, 
microcavity resonances can also be used to create nar- 
row band thermal radiation. [iCl llll | Unfortunately, most 
of the recent designs |8lll2| for perfect absorbers/emitters 
only work for one incident angle and one polarization. So, 
there is a need for wide-angle perfect absorber/emitter 
nanostructures. In this paper, we suggest a structure 
which exhibits a large absorption in the THz regime for 
a wide range of angles with respect to the surface. We 
show that the absorption characteristics are maintained 
even if the uncertainties in the estimated changes in the 
material parameters, due to high temperatures, are con- 
sidered. The proposed structure can be easily manufac- 
tured with today's planar micro-fabrication techniques. 
We also comment on the impact of deviations in the ge- 
ometrical parameters caused by fabricational tolerances. 
The small size of the structure, in comparison to the 
wavelength together with the relatively straightforward 
fabrication, allows for easy integration into various de- 
vices, such as perfect thermal emitters, perfect absorbers, 
bolometers, and very effective light extraction LEDs. 



The suggested structure is shown in Fig.[TJ It consists 
of a metal back-plate (black) with a thickness larger than 
200nm. This is much larger than the typical skin depth 
in the THz regime and avoids transmission through the 
structure. In this case the reflection is the only factor 
limiting the absorption. The thickness of the back-plate 
can be adjusted to the specific needs of the final appli- 
cation, e. g., to obtain good heat transport to sensors or 
to obtain a better stability. On top of the metal plate 
a spacer layer of Silicon-Nitride SiN is deposited with a 
thickness D4 . The structure is terminated by an array of 
metallic stripes with a rectangular cross-section. Their 
arrangement is described by a lattice constant, a, and 
their shape is given by a width, Wuj, and a thickness, 
Wj. In this setup a strong resonance with a large field 
enhancement in the dielectric spacer layer and in between 
the stripes can be obtained, as will be shown later. Ad- 
justing the size of the metal stripes on the top, the cou- 
pling to this resonance can be tuned and the reflection 
can be minimized. 




FIG. 1: Definition of the structure and parameters. An in- 
coming plane wave with an angle 9 to the surface normal 
is considered. The structure consists of an array (lattice con- 
stant a — 2^m) of Tungsten wires (Wt = 0.2/im, W™ = 0.4pm) 
on top of a Silicon-Nitride substrate (Dt = 0.65/xm). The 
Tungsten plate at the bottom must be thicker than 200nm. 
For material parameters see text. 

Due to the scalability of Maxwell's Equations, in prin- 
ciple, the structure can be simulated using dimensionless 



units by dividing all sizes by the lattice constant and us- 
ing oj' = a/A as frequency. However, the Drude model 
used to describe the metal requires frequencies in THz 
and therefore the lattice constant must be assigned in 
the simulation. If a shift in the frequencies of the spec- 
tral features by adjusting the lattice constant is intended, 
a new simulation must be done since changes in the di- 
electric constant would not be considered. 

In the simulation frequency-dependent material pa- 
rameters are required. We calculate them using standard 
methods and adjust their values to take into considera- 
tion the high temperatures. The Tungsten parts (plate 
and stripes) are described by a Drude model 
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with the plasma frequency LOp = 1448TIIz and the col- 
lision frequency oJc — 13THz at room temperature. [13| 
In order to account for the higher temperatures, we use 
an increased value for the collision frequency. Since the 
resistivity of Tungsten increases linearly with tempera- 
ture, a linear dependence of the collision frequency is as- 
sumed, leading to an estimated increase of the collision 
frequency by a factor of 3-5 (wc — 50THz). The plasma 
frequency is assumed to be constant. The simulations 
are repeated for several values of the collision frequency 
(13, 50, lOOTHz) to ensure the stability of the emission 
spectra in a wide range. 

The index of SiN depends on the actual fabrication 
process [ij, |l5|, |l6| and shows a strong wavelength de- 
pendence below 1/Ltm, but is roughly constant for longer 
wavelengths until about 8/im..[17i| For our simulation we 
use experimental values for room temperature obtained 
from Sandia National Lab[l8| and fitted a Cauchy model 
for dielectrics in the infrared to calculate the required 
values. In order to account for temperature increase, we 
assume a linear temperature dependence as found, e.g., 
in Silicafl^l and add a constant value (-1-0.1) to the imag- 
inary part of the experimental data. Again, we perform 
several simulations with different values of the imaginary 
parts of SiN (0.05, 0.1, 0.15, 0.2) for comparison. 

For the simulations we use our own implementation 
of a Fourier-Modal method with a scattering matrix ap- 
proach, also known as Rigorous-Coupled Wave Anal- 
ysis (RCWA).[20|, |2l| Special care is taken in use of 



the correct Fourier-factorization rules to ensure a fast 
convergence. [22| This approach assumes incoming plane 
waves onto a periodic surface defined by the dimension- 
less frequency u)' = uja/27Tc — a/A and the angle 9 to the 
surface normal (Fig.[T]). For the material parameters the 
frequency is calculated using a fixed lattice constant of 
a = 2^TO. The absorption is determined by subtracting 
the total reflection of all diffraction orders given by the 
individual Poynting vectors. The transmission is on the 
order of the numerical error. By Kirchoff's Law the ab- 
sorption equals emissivity and the emission can be calcu- 



lated by multiplying the absorption with the black-body 
radiation, even in the case of structured materials. [3j 

To confirm the results and to study the field dis- 
tributions/energy flow and resistive heating, we use 
the commercially-available finite-element-based program, 
Multiphysics, by Comsol with periodic boundary condi- 
tions in the direction parallel to the surface. In this case 
the absorption is determined by A = 1 — 1 5*211^ — IS'up. 
This program is also used to determine possible eigen- 
modes and their corresponding complex cigcnfrequcncies. 
In this simulation, the material parameters are set to the 
values at the resonance frequency. 

For wide-angle absorption in a small frequency win- 
dow, the center frequency of the absorption peak should 
not shift for oblique incoming plane waves at different 
angles. It was shown before that Tungsten microcavities 
can exhibit a weak angular dependence, 2J, |25[ however, 
the angular range and maximum absorption both are 
smaller than in our case. The additional dielectric layer 
offers the possibility to optimize the cavity geometry and 
the coupling to it lead to a substantial improvement. We 
find that the most important criteria for obtaining a very 
weak angular dependence is to avoid coupling to prop- 
agating — 1*^' diffraction orders. In the simulations we 
choose a lattice constant of a = 2fim, so that the first neg- 
ative diffraction order at large angles {9 « 90°) is given 
by A/a « 2 corresponding to 75TIIz. Any absorption 
peak at higher frequencies couples to several reflected 
diffraction orders and shows a strong angular dispersion 
at the corresponding frequency/ angle combination. For 
a structure with the parameters given in Fig.[Tl we find 
a strong absorption peak for perpendicular incidence ful- 
filling the above criteria. The corresponding reflection 
and absorption spectra are shown in Fig. [51 
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FIG. 2: (Color online) Absorption and reflection for perpen- 
dicular incidence. The Tungsten back-plate is thicker than the 
skin depth and does not allow for transmission. The absorp- 
tion reaches more than 99.9% for u = 69.24THz (A = 4, 3fim). 

The obtained absorption peak has a center frequency 
of i'o = 69.24THz with a half-width at half-maximum of 
10.27TIIz corresponding to 14.8% of the center frequency. 



Simulations with Multiphysics showed, the field is able 
to couple into the dielectric layer in contrast to regions 
outside the resonance. A strong enhancement of the elec- 
tric field can be found in the region between the metallic 
stripes mostly located in the dielectric (Fig. [3]). Numeri- 
cally, we find an eigenmode showing the same field pat- 
tern at an eigenfrequency of 70.05THz (real part of the 
complex eigenfrequency). In this simulation the mate- 
rial parameters are fixed to the values at lyQ. Changing 
the thickness of the metal plate or the air space in front 
of the structure changes the eigenfrequency by less than 
O.OlTHz, once the air part in front of the structure is 
chosen large enough. Although losses in the dielectric 
layer occur due to the imaginary part, the energy flow 
described by the Poynting vector reveals that the main 
absorption takes place at the surface of metallic back- 
plate and the rear side of the metal stripes. 
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FIG. 3: (Color online) Magnitude of the z-component of the 
electric field at the resonance. The linear scale ranges ranges 
from 0.0 (light green) to 25000 (dark red). The maximum of 
the field is located between the metal stripes in the dielectric. 
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FIG. 4: (Color online) (a) Angular peak dispersion of the ab- 
sorption peak. The additional peaks around 110 and ISOTHz 
show a strong dispersion although they are caused by cavity 
modes as well, (b) For to 40° the absorption is above 99.7%. 
Even at 70° (80°) still 83.3% (58.1%) of the incoming energy 
is absorbed. 



If the structure is used as a perfect absorber, e.g., in 
sensor applications, it is important to absorb as much 
of the radiation as possible, independent of the direction 
of incidence. In Fig. 3] we show the angular dispersion of 
the peak together with the spectra for several angles. For 
angles up to 40° the peak becomes slightly narrower and 
the center frequency increases by O.STHz, but the max- 
imum absorption is above 99.7%. For larger angles the 
peak starts to drop, but is still above 80%, even for waves 
impinging on the structure with an angle oi 9 — 70°. In 
Fig.mja), two additional absorption peaks at higher fre- 
quencies are visible. Both peaks are caused by resonant 
modes in the dielectric as well, but their angular depen- 
dence shows a very strong dispersion. If the parameters 
are chosen so that the lowest peak has a higher center 
frequency (above 75TIIz for the given lattice constant), 
it shows a weak angular dispersion until the — l'^' starts 
propagating and then the absorption vanishes quickly. 

Finally, we compare the emission properties of our 
structure with a regular black-body emitter in Fig. [51 We 
plot the black-body emission using Planck's law at a tem- 



perature of 1176K, so the maximum of the emission cor- 
responds to the center frequency of the absorption peak. 
The emission of the structure given by the product of the 
emissivity and the black-body radiation is also plotted in 
Fig. [5] as a solid line. The combined effect of the two as- 
pects suppresses the sides of the peak even stronger, so 
a selective emission is achieved. 

For these high temperatures, the material parameters 
are not known exactly and must be estimated. Since 
estimates are always a possible source of errors, we en- 
sured that the spectral properties exist in a wide range 
of parameters. We tested the results for different ad- 
justments in the collision frequency of Tungsten {coc =13, 
50, lOOTHz) and the imaginary part of the SiN dielectric 
(0.05, 0.1, 0.15, 0.2). For all combinations of these pa- 
rameters, the absorption reaches values of 93% or more, 
except for the case with ujc — 13TIIz and Ime = 0.05 for 
which only 82% are achieved. The center frequency of 
the absorption peak remains unchanged. However, an 
increase in the imaginary part of SiN leads to an broad- 
ening of it. In general, an increase in the collision fre- 
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FIG. 5: (Color online) Black-body radiation (dotted), ab- 
sorption (dashed), and emission (solid) of the structure, if 
the temperature is adjusted so the maxima of the black-body 
spectra and the absorption peak coincide using Wien's law 
(T=1176K). Emission and absorption are roughly the same, 
due to high value of the black-body radiation. 



quency allows for a wider range of possible adjustments 
in the imaginary part of the SiN. 

Although we do not present the data in this paper, we 
also conducted several simulations to determine the sen- 
sitivity of the absorption peak to fabricational tolerances. 
In this study, we varied all 3 parameters (W^,, Wt, Df) 
by ±50nm individually and combined. In all cases an ab- 
sorption of more than 98,8% was obtained. Wider (nar- 
rower) stripes lead to a small increase (decrease) in the 
frequency of less than ITHz. For thicker stripes, the cou- 
pling became weaker and the absorption dropped by 1%. 
Changes in the thickness of the dielectric layer caused 
the strongest change in the absorption spectra, shifting 
the center frequency of the peak by about ±4THz. Since 
the field was concentrated mainly in this area, deviations 
in the spacer thickness changed the effective size of the 
cavity and, therefore, the frequency of the cavity mode. 
This makes it important to control the thickness of the 
dielectric layer very precisely in the fabrication process. 

In summary, we presented a design for a perfect ab- 
sorber/thermal emitter with a very high absorption peak 
over a wide range of angles. We demonstrated that the 
absorption peak is stable with respect to the estimated 
changes in the material parameters at high temperatures. 
We also studied the impact of deviations in the geome- 
try caused by fabricational tolerances. From this it can 
be expected that the spectral features are also present in 
realistic samples. Since the structure can be built with 
planar fabrication techniques, it offers an interesting ap- 
proach to wide angle perfect absorbers/emitters. First 
studies using a lattice constant of several mm showed 
that it is possible to obtain similar effects also for fre- 
quencies in the GHz range. Due to the wide angular ab- 
sorption this can offer a new way to avoid reflections in 
microwave experiments with an easy-to-build structure 
on length-scales smaller than the wavelength. 
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